Run plot_mech_bands.m